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Abstract 

Employing a variant of the Multicanonical Ensemble we study 
twisted spin configurations on periodic boxes in the D = 2 0(3) non- 
linear sigma model for /3-values inbetween 1.55 to 3.1. The free en- 
ergy difference of twisted spin configurations is determined from the 
constraint effective potential. The finite size scaling behavior is in ac- 
cordance with the asymptotically free nature of the continuum theory. 
Upon certain reasonable assumptions we determine the A/3(/3)-shift of 
the stiffness correlation length £ s . The mass-gap as determined by our 
analysis is rriQ = 79.6(1.9) A; #. This value agrees with the analytical 
result of the thermodynamic Bethe Ansatz tuq = 80.1 A/ att . 



1 Introduction 



The two dimensional 0(3) nonlinear sigma model is a asymptotically free 
field theory pQ. It deserves special interest in view of asymptotically free 
four dimensional nonabelian gauge theories describing the strong interactions 
inbetween Quarks. A fundamental feature of these theories is the existence 
of a nonvanishing mass-gap rriQ. Recently the mass-gap of the sigma model 
has been calculated analytically on the basis of the thermodynamic Bethe 
Ansatz E] . Choosing the simplest regulator, namely a two dimensional 
hypercubic lattice with lattice spacing a, and the standard nearest neighbor 
action functional with nearest neighbor coupling f3 we consider 

s = -p £ , (i) 

<i,j> 

where the fields (j)f are 3-component unit vectors and the sum < i,j > runs 
over the nearest neighbors on the lattice. The mass-gap correlation length 
£o = ttiq 1 has the value 

I-tt/2 2tt/3 qoi 

corresponding to the exact mass gap m = 80.0864 Ai aU [H], Ai att denoting 
the 3-loop lattice cut-off parameter [1]. The past numerical studies of the 
path integral Ej consistently overestimated this number. These studies 
considered infinite volume mass-gap correlation length values £ U P to about 
hundred lattice spacings a with the help of cluster algorithms. It was also 
attempted to enlarge the accessible correlation length region utilizing real 
space renormalization group methods 5J. It was found, that asymptotic 
perturbative scaling only sets in at mass-gap correlation length values, which 
can hardly be fitted into the memory of todays computers. In addition it 
may also be appropriate to mention the still standing criticism of Patrascioiu, 
Seiler and others on the continuum limit of the sigma model 0. These 
authors argue in favor of a the existence of a phase with vanishing mass-gap 
in the sigma model, which is however not supported by the numerical data 
jH]. It is therefore a challenge to confront the known analytical result on 
the mass-gap with the precision numerical simulations of the path integral. 
From such a study we might also gain better insight on how to accurately 
determine the mass-gap of lattice QCD. 
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Finite size scaling theory 9\ has been very successful throughout the 
years in predicting the behavior of the infinite volume correlation length of 
spin systems at criticality from the properties of finite systems, which in 
this paper are taken to be square boxes with linear extent L and periodic 
boundary conditions. In a situation where the control parameter xq = ^q/L 
is large and, even at the critical point, where the correlation length is infinite, 
the path integrals singular part of the free energy on finite volume systems 
is a function of the control parameter x alone. The Fisher scaling analysis 
leads to reliable determinations of the correlation length divergence critical 
exponent v in statistical mechanics models. 

In this paper we attempt a generalization of this idea to the D = 2 
nonlinear sigma model. Its asymptotic freedom fixed point is located at 
infinite value of the bare nearest neighbor coupling. The quantity under 
consideration is a free energy difference of twisted relative to untwisted spin 
configurations - the spin stiffness - forming a Bloch wall, whose finite size 
scaling behavior is dominated by logarithmic terms in the control parameter 
Xs = £,s/L in accord with the perturbative treatment of the path integral 
in the one-loop approximation. The quantity £ s hereby denotes the stiffness 
correlation length. We will argue, that a precise extraction of the A/3(/3)-shift 
corresponding to the stiffness correlation length is feasible up to j3- values as 
large as 3.1, upon reasonable assumptions on the coefficients of the finite size 
scaling law. The integration of the shift in junction with a start value for 
the mass-gap correlation length £o allows the determination of the mass-gap 
m . The paper is organized as follows: Section 2 is devoted to the theoretical 
background and introduces the constraint effective potential, the considered 
free energy difference and discusses the finite size scaling hypothesis. Section 
3 contains details of the Multicanonical Ensemble simulation. In section 4 
we present our numerical analysis and results. Section 5 concludes the of the 
paper. 

2 Theoretical Considerations 

In statistical field theories, which share the universality class of continuous 
ferromagnets, it is well known on the basis of the e-expansion and the renor- 
malization group JU] , thai the critical behavior of the system can be studied 
in the symmetry broken phase of the theory by considering path integrals 
of configurations, which interpolate inbetween regions of different order pa- 
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rameter orientation. The helicity modulus Y carries the information on the 
nonanalytical behavior in the symmetry broken phase of the theory. It is 
common to define the helicity modulus by the response of the system with 
respect to a twist angle 0, which is applied microscopically to the fields on the 
boundary of the lattice in one of the lattice directions. Denoting F = —InZ 
the free energy, the free energy difference 

AF Y {0) = F{8) - F{9 = 0) (3) 

then defines the helicity modulus [U] 

which is finite in the thermodynamic limit and allows for a geometry in- 
dependent i.e., O-independent, characterization of the Bloch wall. Hereby 
denotes Q the cross section orthogonal to the direction of the twist. E. g., 
on a hypercubic lattices with linear extent L one has Q oc L D ~ l and one 
expects AFy(0) oc L D ~ 2 . The critical behavior of the theory in terms of 
the helicity modulus is expressed by Josephsons scaling law [12 : Y = Ayt^ J 
with fij = (D — 2)u. It is clear, that the presented scaling laws refer to the 
ferromagnetic case. In fact, when the dimension D = 2 is approached from 
above for field theories with continuous symmetry , like the O(N) nonlinear 
sigma model with iV > 3 0, the ferromagnetic character of the theory is lost 
and the helicity modulus Y according to eq.Q vanishes for any value of the 
bare coupling (3. In dimension D = 2 the infinite volume system exhibits do- 
mains of different order parameter orientations, which are free to fluctuate. 
This is the content of the Mermin Wagner theorem. 

At the heart of the problem is a study of the correspondingly defined free 
energy difference AF s (0) on finite lattices in the D = 2 0(3) nonlinear sigma 
model at large values of the nearest neighbor coupling i.e., in a situation in 
which the mass-gap correlation length £o exceeds the linear extent L of the 
symmetric hypercubic lattice. In this situation one expects AF S (9) to be 
controlled by x s = £ S /L, where £ s denotes the stiffness correlation length. 
Based on the renormalization group, the e = D — 2-expansion and the 1- 
loop approximation, the calculation has been performed for the case of a 
continuum square box with fixed twisted boundary conditions [THIEI]. We 

1 We refrain from a discussion of the 0{2) symmetric XY model. 
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quote the result 



AF S (Q) = —[Infa/L) + Mnfo/L)] + R(Q) + C(oo). (5) 

A few remarks are in order: 

1) According to the expected Fisher scaling AF S (Q) splits into a singular 
part, which is a function of x s = £ s / 1 L alone and a regular contribution R{&). 
The stiffness correlation length £ s follows the perturbative renormalization 
group, and thus in the large /3-limit has the form 

UP) oc (6) 

However, £ s is not identical to the mass-gap correlation length £o as is the 
control parameter x s as compared to Xq. The stiffness correlation length £ s 
describes the crossover of the system from its small volume behavior at large 
values of x s into a state of various domains of order parameter orientation 
on large volumes. It is expected to be larger than the mass-gap correlation 
length £o- Recently Billoire performed a high statistics numerical calculation 
of the stiffness correlation length at /3-values far in the perturbative 
region f3 = 5, 10 and 20 respectively. Expressing the measurement in units 
of the exact mass-gap correlation length we quote £ s = 9.39(1), 9.46(1) and 
9.47(l)£o heor - In this paper we assume £ s takes its perturbative value that at 
a /3-value of 3. In addition we expect a region of /5-values, where £ s as well 
as £ are governed by the same beta-shift A/3(f3). Let us denote with A/3(/3) 
the change of the coupling (3 to (3 — /3 — A/3(/3) corresponding to the decrease 
of the stiffness correlation length by a factor of 2 : £ s (/3) = |£ s (/3). 

2) The free energy difference AF S (G) at = defines the G independent 
spin stiffness p = AF S (6). The large (3 perturbative limit of p s = ~q§j;P 
defines the spin stiffness constant whose value is Precision numerical sim- 
ulations of the spin stiffness ^Hl QHj confirm this theoretical expectation. In 
the context of the present paper we assume that the spin stiffness coefficient 
proportional to InL takes its perturbative value in the whole considered (3 
region. 

3) Upon insertion of eq.fJOJ) into eq.© we observe, that the leading term in 
AF S is j3Q 2 /2. This is the classical action difference of a twisted configuration 
at twist angle G relative to a configuration with twist angle = 0. The 
dominant effects of quantum fluctuations at fixed (3 are linear in InL and are 
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given by a term —Q 2 lnL/A7T. They drive the spin stiffness to smaller values 
on the larger system sizes, as expected. 

4) The regular contribution R(Q) in eq.(0) can be calculated in the small 
expansion. For the boundary conditions considered in one finds 
R(Q) = —2.501^4 + (9(0 6 ) i.e., a negative contribution at finite 9 values. 
We note that these contributions are nonuniversal in character and depend 
on the details of the implementation of the twist. 

Numerical simulations of free energies are computationally difficult. In 
the standard approach one integrates the expectation value of the action 
difference of systems with nonzero twist G and twist = along the nearest 
neighbor coupling parameter direction (3. One can avoid the integration of 
the action by differentiating the path integral with respect to at a value 
of = P^ ll7j. In this paper we consider the constraint effective potential 
(CEP) [TH] of the mean field of the theory. 

The consideration of the CEP is motivated by the analogous problem aris- 
ing in the interfacial case in Z(2) symmetric theories and recently numerical 
simulations of it have been conducted with the help of Multicanonical En- 
semble simulations [TH] . One considers the mean field M = jp J2i Si of the 
Ising fields. Its probability distribution P(M) contains at the same time 
the information concerning bulk as well as interfacial properties of the sys- 
tem. Relevant for interfacial effects are states at M = on finite boxes 
with periodic boundary conditions. In this region of the phase space config- 
urations contain two domains of opposite order parameter orientation and 
two interfaces are formed. Widoms scaling law can be analyzed [2U|. Lee 
and Kosterlitz have performed a finite size scaling analysis of the CEP at 
criticality [21] leading to a determination of v. 

In case of the theory with continuous 0(3) symmetry we introduce the 
3- component field 



3 a = 7iFX>; 



(7) 



and its absolute value, denoted the mean field $ in the following 




(8) 



which will have the probability distribution 



P($) oc ¥e~ u ^ 



(9) 
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in the canonical ensemble of the theory. The function U(<&) is the CEP of 
the theory. It can be obtained by rewriting the partition function 

Z = J D0e /, £»."* Sf * s +» (10) 

into r°° 

Z = / d$&e- u <*\ (11) 
J o 

which can be achieved by introducing suitable (^-functions and integration 
upon the remaining degrees of freedom. We note that U(3>) is defined up to 
constant, which can easily be absorbed into a multiplicative normalization of 
Z. In addition a factor $ Ar_1 at iV = 3 appears in the path integral definition 
of U. This factor is proportional to the surface of a iV-sphere of radius $ 
and accounts for the degeneracy of states with respect to O(N) rotations. 
Without taking this phase factor properly into account the CEP is a singular 
function at $ = 0. 

The CEP attracted in the past attention in the context of the Higgs- 
models represented by O(N) symmetric theories in D = 3 and D = 4. There 
it was studied with analytic methods [22] as well with numerical simulations 
|23j . In these models the CEP has at a finite value 5> m j„ a minimum, which 
corresponds to the Higgs field expectation value. We may note, that both an- 
alytic, as well as the numerical considerations in these models were concerned 
with the shape of the CEP in the vicinity of its minimum. 

In the context of the present paper we are however concerned with the 
CEP at values of the mean field <£> equal to $ = 0. We argue that states at 
$ = correspond on hypercubic boxes with periodic boundary conditions to 
Bloch walls carrying formally a twist angle of = 2ir. This can be easily 
seen on the classical level by noting, that the field configuration 

1 = cos (—±) 4> 2 = sini— 1 ) 4> 3 = (12) 

Lj Lj 

minimizes the action functional for the given periodic boundary conditions at 
$ = 0. This configuration exhibits on the scale of one lattice spacing a spin 
twists of value ^ in the 1-direction of the lattice, which upon integration 
in the 1-direction adds up to a total twist angle of 6 = 2n. We expect this 
configuration and the added quantum fluctuations to dominate the $ = 
state. Thus in the framework of the CEP it is strongly suggested, that the 
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singular part of AF S (Q = 2tc) exhibits the same finite size scaling as the 
constraint effective potential difference, or potential barrier 

AU = U($ = 0) - U($ min ) = AF S (Q = 2tt) + R, (13) 

up to regular terms R, provided the CEP exhibits a maximum and mini- 
mum corresponding to the twisted and untwisted states. It is this finite size 
scaling hypothesis, which we will examine in the subsequent analysis of this 
paper. We will find, that the numerical evaluation of the CEP confirms the 
hypothesis. 



3 The Multicanonical Ensemble Simulation 

By the use of standard Monte Carlo algorithms it would be very difficult 
to sample values of $ close to zero, if the value of the nearest neighbor 
coupling f3 is larger than about 1.6 on lattices of reasonable linear size L. To 
overcome this difficulty, we modify the importance sampling by introducing 
a Multicanonical- weight factor into the Monte Carlo sampling process [TH] . 
We remark, that the idea of Multicanonical Ensemble simulations consists in 
modifying the Boltzmann-weight for the purpose of the improvement of the 
actual Monte Carlo sampling process. The modification is however done in 
such a way, that the effect of the Multicanonical-weight will be removed in 
a well controlled way. Thus in the end the CEP can be determined in the 
canonical ensemble of the theory. 

In case of the CEP the Multicanonical-weight factor is chosen to be a 
function of the mean field 5> evaluated on each single configuration and will 
be denoted W mc ($). It is in principal defined on the real interval [0, L 2 ] 
and it is sensible to split the Multicanonical-weight into a part, which takes 
care of the degeneracy with respect to 0(3) rotations, and a yet unknown 
contribution W mc (3>). Thus for the purpose of the Monte Carlo simulation 
we consider the Multicanonical-weight 

W mc (l>) = -2 ln$ + W mc ($) (14) 

and 

e -s+w mc (5) ^ 15 ) 

to be the Boltzmann-factor, which generates the Markov process. In order 
to obtain a representation of the Multicanonical-weight W mc , which in a 
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computer is easily calculable, we choose W mc to be a polygon, which on 
a finite set of i = l,m intervals Ii : < $ < for each interval is 

characterized by two parameters gi and hf. 



One may view the parameters hi as magnetic sources, which are applied on 
piecewise intervals of the operator $ to the theory. In one of our earlier 
publications we therefore named the resulting ensemble the Multimagnetical 
Ensemble |20j . We found it sufficient in the actual numerical simulations, 
to first determine a value $ max in a standard simulation. In this situation 
Qmax is calculated in such a way, that all the relevant structure of the CEP, 
namely the location of its minimum, and the states at $ = are included in 
the considered 5>-interval. It means that we will only evaluate the shape of 
the CEP for values < $ < <& max - We then choose for reasons of simplicity 
an aequidistant partition of the considered ^-interval into m = 25 or m = 40 
different intervals Ii, which in the actual simulations turned out to be a 
sufficiently fine polygon approximation to the Multicanonical-weight W mc . 

Multicanonical Ensemble simulations are defined by the requirement, that 
the considered operator exhibits an almost constant, ideally a constant, prob- 
ability distribution in the simulation. If in the canonical ensemble the den- 
sity of states function with mean field $ is denoted n(i>), then the ideal 
Multicanonical-weight factor is 



Though trivial, we note, that a priori the Multicanonical-weight is not known. 
Less trivial we remark, that it is possible to obtain an estimate of the weight 
factor from the Monte Carlo simulation. In the vicinity of a given value of 
the operator $ one may sample the phase space in an attempt to estimate the 
density of states. In practice we perform for each of the considered intervals 
Ii a separate simulation with 5>-values constrained to the given interval. In 
the simulation each parameter hi is then determined under the requirement 
that the probability distribution of the mean field $ approximates the Mul- 
ticanonical distribution at its best. The full set of all parameters hi serves 
the purpose of the construction of the complete Multicanonical-weight ac- 
cording to eq.(JT3J) and eq. lfTo'j) . Clearly, for given hi the parameters Qi can 
be chosen in such a way, that the resulting polygon is a steady function of 
$. We note, that one parameter value, e.g. g± is left undetermined by this 



W mc (<$>) = gi + hi®. 



(16) 




(17) 
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procedure. This again corresponds to the overall normalization freedom of 
the path integral. In our actual simulations we have found, that this sim- 
ple and robust procedure results into acceptable Multicanonical-weights and 
distribution functions. 

Once the Multicanonical-weight factor has been determined, we imple- 
ment one sweep of the Monte Carlo sampling by a combination of a Swendsen- 
Wang reflection cluster update |2HES1, followed by a subsequent accept-reject 
decision, which then depends on W mc alone. It means, that we define the 
Cluster degrees of freedom from the cinetic term of the action, and, that 
we replace the usual equal probability rule of the Cluster update, by a 4- 
hit Metropolis accept-reject decision in all the Cluster degrees of freedom, 
depending on the magnetic properties of the system. At this value for the 
number of hits we obtained average acceptance rates for moves of the system 
at about one half for almost all the considered (3 and L values. In order to 
monitor the performance of this update procedure, we divided the considered 
^-interval into an arbitrarily chosen 8, but equally spaced bins. We measure 
the average flip-autocorrelation time 77 for a transition of the system from 
the first to the eights bin, or vice versa. While due to the partitioning some 
arbitrariness in the definition of tj is induced, we nevertheless expect, that 
it will exhibit the basic aspects of the Monte Carlo time dynamics. Fig. 1) 
contains data for the quantity 7/ obtained on a L = 36 lattice for /5-values 
in between 1.6 and 3.0 (circles), and for a fixed /5-value (3 = 2.4 on lat- 
tice sizes ranging in between L = 20 to L = 68 (triangles). One notices 
a typical time scale of several thousand sweeps, which in general is char- 
acteristic for our simulation. In detail we observe, that at a given /5-value 
the flip-autocorrelation time stays almost constant, or even decreases with 
increasing lattice sizes. We attribute this favorable property to the nonlocal 
nature the algorithm. We also observe on the given lattice a rapid increase of 
the flip-autocorrelation time with increasing (3. Such an increase is expected, 
as the simulational complexity of the theory will rise in the proximity of the 
asymptotic freedom fixed point. 

In the Multicanonical simulation we measure the Multicanonical proba- 
bility distribution P mc ($) of the mean field $. It is related to the probability 
distribution in the canonical ensemble P($>) by a simple reweighing step 



and thus with the help of eq.© the CEP can be determined in the Multi- 
canonical simulations. Throughout our simulations we have fixed the mini- 




(18) 
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mum value of the CEP to be zero. In addition we perform a bias corrected 
jackknife error analysis for all the measured quantities. The typical statistics 
accumulated for each data point in our simulation was inbetween 1 and 5 
megasweeps, depending on the considered /5-value. 

4 Numerical Analysis and Results 

The mass-gap correlation length £ of the D = 2 0(3) nonlinear sigma model 
has been calculated in previous studies at not too large values of (3 and we 
quote here the results of Wolff 0: At the /3-values (3 = 1.6,1.7,1.8 and 
1.9 i Q /a takes the values 19.07(06), J 34.57(07), 64.78(15) and 121.2(6) respec- 
tively. In this paper we study small systems of sizes ranging in between 
L = 8 up to L = 82. For /3-values larger than (3 = 1.6 the mass-gap 
correlation length £ is comparable with the considered lattice sizes. Ac- 
cordingly it is expected that the stiffness correlation length £ s exceeds the 
linear size of the systems and the onset Fisher scaling may be expected. 
In a first set of simulations we have studied the L-dependence of the CEP 
at twelve values of the nearest neighbor coupling, namely at the /3-values 
13 = 1.6, 1.7, 1.8, 1.9, 2.0, 2.1, 2.2, 2.3, 2.4 and at the (3- values (3 = 2.6, 2.8 and 
(3 = 3.0. In a second set of simulations the /9-dependence of the CEP po- 
tential was determined on a L = 36 lattice at /3-values ranging inbetween 
f3 = 1.55 up to (3 = 3.1. In total we have accumulated 468 pairs of couplings 
f3 and lattice sizes L within our simulation. 

It is important to check our basic assumption, that states at mean field 
$ = carry a twist angle of = 2tt. For this purpose we use a special 
graphical representation of the field configuration, which is exhibited in Fig. 
2a) and Fig. 2b). On the given configuration we apply in a first step a global 
rotation which equals the arbitrarily chosen 3 component of the field $ a to 
zero. We then map the 1 and 2 components of the fields along paths in the 
lattice onto the complex plane. In this mapping we connect fields, which are 
neighbors within the path, by a line. Fig. 2a) contains all such mappings, 
which are obtained by considering all L paths along the main axis of the 
lattice in the 1-direction. Fig. 2b) contains the analogous graphs for paths 
along the main axis in the 2-direction. The considered lattice size is L = 36 
and the nearest neighbor coupling is (3 = 8, far in the perturbative region. 
The value of $ for the considered configuration was $ = 0.003904, which 
is very close to zero. One clearly observes the winding angle of 27r for the 
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latter case i.e., for all paths on the main axis of the lattice in the 2-direction, 
which is clear evidence in favor of the above assumption. It is noteworthy to 
remark that fluctuations around the classical state eq.(fl2^1 are sizable even 
for such a large /3- value. They increase as (3 is lowered and decrease as (3 is 
enlarged. In addition we also mention, that states at $ = apparently break 
the cubic group lattice symmetry i.e., in the here considered configuration 
the twist appears in the 2-direction of the lattice. We have checked, that 
statistical independent configurations at $ = assume the possible different 
twist directions with equal probability. From a theoretical point of view it 
must be expected that this additional degree of freedom, as compared to a 
theory with a fixed direction of the twist, contributes to the unknown regular 
contributions of the finite size scaling relation eq. QlHj) . 

In Fig. 3a) and Fig. 3b) we display the CEP, as obtained by our simula- 
tions on lattice sizes L = 24, 36 and L = 70 and at a values of (3 — 1.6 and 
(3 = 2.4. 

At (3 = 1.6 we observe a crossover from a situation, where for the smaller 
L- values the CEP for states at $ = exhibits a barrier AU > 0, to a situation 
at the largest L-value, where the barrier vanishes. Thus the decrease of the 
spin stiffness with increasing L is observed. This is attributed to the quantum 
fluctuations around the classical state with twist angle G = 2tt. We may 
remark, that in simulations of ferromagnets, namely the D = 3 0(3) sigma 
model in its symmetry broken phase, we have witnessed an increase of the 
potential barrier with increasing L. 

At 13 = 2.4 the mass-gap correlation length and £ s exceed the linear lattice 
size by a large factor. Correspondingly a crossover cannot be observed and 
AU ~ 20 is finite on the considered range of lattice sizes. Again its decrease 
with increasing L is witnessed. The large value of AU corresponds to an 
exponentially large suppression of states with mean field $ = 0, even if the 
phase space factor of eq. © is divided out from the probability distributions. 
In fact, a simulation of the D = 2 0(3) nonlinear sigma model at (3 = 2.4 on 
small volumes resembles, as far as order parameter orientation is concerned, 
on the first sight the behavior of a ferromagnetic system. It is only, that 
on lattices with linear extent comparable and larger to £ s , the system will 
approach its true vacuum state. 

From the shape of the CEP we have carefully determined the quantity 
AU = [/($ = 0) — U(<& min ). For this purpose we employ fits to the shape 
of the CEP in the vicinity of the minimum at $ m i n and in the vicinity of 
$ = 0. In the vicinity of $ = and for the determination of U(3> = 0) we 
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describe the CEP by a parabola. Such a functional form can be expected on 
the basis of a classical argument 



U = [/($ = 0) + «i$ 2 (19) 

and is consistent with the data. We use it in a fit range of $-values close 
to $ = 0, for which corresponding Xdof~ values of the fits are smaller than 
unity. Inspecting the shape of the CEP in vicinity of its minimum we note 
an apparent asymmetric behavior and therefore we use the form 

U = U($ min ) + - <J mm ) 2 + &($ - $ mm ) 3 (20) 

for its analytic description and the determination of U($ m i n ). Analogous 
remarks as in the previous case on the fit-intervals and the xl. ./.-values 
apply. 

We expect, that discretization effects of the lattice theory, as compared 
to the continuum, lead to correction terms to the spin stiffness. The leading 
contributions are of the order I /a 2 . On the lattice and in the classical ap- 
proximation we calculate the action difference ASojatt of a field configuration 
with twist = 2-7T relative to untwisted fields. It has the expansion 



AS 0Mt = AS ML) + 0(- 4 )}, (21) 

a 



with p{L) given by 



7T 2 



M = 1 - C 'W4>- (22) 

Classically the constant C g is equal to unity, and AS , C = 2ir 2 (3 denotes 
the action difference of the Bloch wall in the continuum. We will assume 
here, that the possible form of the 1 / a 2 -corrections is already determined on 
the classical level. Correspondingly we assume in the fully fluctuating theory, 
that additional contributions can be accounted for by a nonunity value of the 
parameter C g . The potential barrier in the continuum AU C is then related 
to AU evaluated on the lattice via the relation 

AU C = P (L)- l AU, (23) 

which then allows the comparison of the lattice data with the continuum 
scaling form. 
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We have analyzed the finite size scaling of almost all of our AU data by 
means of one x 2_ fit with the form 



AU=[1- C g ^—^][irln^ 8 /L) + Alnln(UL) + R]. (24) 

Aiming at a precise determination of the mass-gap £o we implement our 
knowledge about the perturbative limits of the spin stiffness p and the spin 
stiffness correlation length £ s . We constrain certain parameters and adopt 
the following scheme: 

1) In accord with theoretical expectations the prefactor of the single log- 
arithmic scaling term oc ln(£ s /L) of the continuum free energy is fixed to its 
exact value ir. This is expected on the basis of our previous discussion on 
the spin twist and the spin stiffness. We leave the prefactor A of the double 
logarithmic scaling term to be a free parameter. It may be tolerated, that 
some of the omitted 0(1) higher order terms of the loop expansion in eq. (|24j) 
may be represented effectively by a value of A, which differs from tt. 

2) The quantity dpln^ a ((3) is expanded in inverse powers of the nearest 
coupling parameter (3. Thus 

d.lnUP) =^-^ + °-^ + (25) 

represents an analytic form, which incorporates the known 3-loop behavior 
of the mass-gap correlation length. It is valid, if in the asymptotic scaling 
region £ and £ s follow the same (3 function. Without loss of generality we 
have included 4 free parameters jk for k = 3, 6. These additional degrees 
of freedom allow for deviations from perturbative scaling in the small /3-value 
region, where the dip in the Af3(f3) function is expected. 

3) Given a start value for the stiffness correlation length at a value of the 
nearest neighbor coupling (3* in the perturbative region of the theory, £ s can 
be integrated to smaller values of (3 via 

InUP) = f dplntsiP) + InUP*), (26) 

if within the fit the quantity dpln^ s (0) is represented by its expansion eq.(|25j). 
We exploit the recent measurements of the spin stiffness correlation length in 
the asymptotic scaling region ^Hj and fix the integration constant ln^ s (/3*) 
by choosing (3* = 3 and £ s (/3*) = 9. A7^ heor (/?*). Note then that at (3* = 3, 
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£ s has a value of 934051a. Any significant deviation of £ s from its measured 
value in the asymptotic scaling region would appear unreasonable to us at 
such large values of the correlation length. 

4) We also express the functional form for the quantity C g (j3) through an 
expansion in inverse powers in j3 

k=i p 

upon introduction of 4 free parameters 5k with k = 1,...,4 parameterizing 
the departure of I /a 2 corrections from the classical result. At infinite value 
of (3 the classical result is recovered. 

5) The onset of finite size scaling is expected for large values of the control 
parameter x s » 1. Into the actual fit we include data from the /3-interval 
1.65 < (3 < 3.1, leaving us with 428 independent data points. For these data 
and for the considered lattice sizes our final fit only includes data, whose 
control parameters obey the inequality x s > 4.7. Such large values of the 
control parameter x s hopefully put us in a region of the theory, where we 
can reliably trust our finite finite size scaling ansatz eq. ()24j) . 

The resulting 10 parameter x 2- ht for the parameters A, R, 7^ and 5k with 
k = 1, 4 has been executed with the double precision version of the CERN 
libraries MINUIT program. The error analysis for the parameter values and 
for all other derived quantities of this paper has been performed by a repeated 
execution of the fit on several data sets. In each data set any datum in the 
fit is distributed gaussian around its mean with a variance corresponding 
to its error. The fit has an excellent Xd.o.f~ Y8 ^ ue °f Xd.o.f = 0-75 and is 
displayed in Fig. 4) and Fig. 5) in comparison with the measured AU values. 
The fit as denoted by the lines in the figures captures the /3-dependence of 
the data in Fig. 4) as well as the L-dependence displayed in Fig. 5). As 
expected by the theoretical reasoning the fit supports a negative regular 
contribution to the scaling law eq.(J2H) R = —3.45(20). It can be compared 
with the mentioned result of the small B-expansion of (9(G 4 ) on boxes with 
fixed boundary conditions at a formal value of the twist angle 2n, which 
is —5.002. However there is no theoretical argument, that both numbers 
should be identical, rather they should differ by contributions induced by the 
different boundary conditions. We also find a positive coefficient A = 1.79(7) 
attached to the double logarithmic scaling terms. The former value differs 
from the expected expected value tc, which we attribute to the subleading 
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nature of the double logarithmic terms. In Fig. 6) the fit result for the 
constant C g is given. Finite grid effects are largest for /^-values around (3 = 
1.9 and exceed their classical value in the constant C g by a factor of about 
3 there. With increasing (3 a very slow approach to the classical result is 
witnessed. In Fig. 7) the logarithm of the stiffness correlation length as 
well as its /3-derivative are displayed. Note again that we have fixed £ s at 
(3* = 3 to the measured perturbative result, the fat circle in the figure. The 
calculated /^-derivative of ln£ s as displayed by the solid curve in the inlay of 
Fig. 7) exhibits a very clear "peak like" deviation from 3-loop asymptotic 
perturbative scaling, the dashed curve in the inlay of Fig. 7). The peak 
corresponds to the dip in the A/3(/3)-shift as noticed in the earlier studies 
of the sigma model. It is located at a /3-value of about (3 ~ 1.75. Notably 
our data demonstrate, that deviations from 3-loop asymptotic scaling are 
present up to (3- values of about (3 = 2.5, while for all larger values of (3 
asymptotic perturbative scaling is realized within the numerical precision. 
This finding is a genuine outcome of the present analysis and also applies 
to the A/3(/3)-shift corresponding to the decrease of the stiffness correlation 
length by a factor 2. The A/3(/3)-shift is displayed in Fig. 8) for /3-values 
larger than (3 = 1.8, the solid curve in the figure. It is compared with the 
3-loop asymptotic perturbative scaling result, the dashed curve in the figure 
and with results from Hasenfratz et. al [Sj. Table 1) collects few entries for 
the shift at selected values of f3. We thus observe a somewhat slower approach 
to asymptotic scaling as it was previously anticipated in the the numerical 
simulations. A similar observation has already been reported earlier [2*o] . 

The calculation of the mass-gap of the theory in units of the 3-loop lattice 
cut-off can be attempted under the nontrivial assumption, that mass-gap as 
well as stiffness correlation length share a common /3-value region, in which 
their flow towards the asymptotically free fixed point with varying couplings 
is governed by a universal A/3(/3)-shift. In this scaling region the ratio £o/£s 
stays constant. It is common believe, that this region stretches beyond the 
region where asymptotic scaling is observed. Providing a start-value of the 
mass-gap correlation length at a small value of f3 outside the perturbative 
scaling region, but within the so called nonasymptotic scaling region, it may 
be feasible to integrate the mass-gap of the theory with the help of d^ s (f3) 
up to a /5-value (3* in the asymptotic scaling region. The existence of a 
such a universal behavior beyond the asymptotic regime is however by no 
means guaranteed. Referring to the results of the preceding paragraph we 
are confident now that the choice (3* = 3 moves us into the asymptotic scaling 
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region. For the start value of the mass-gap correlation length we choose one 
of Wolffs numbers 

£ (/3 = 1.8) = 6478(15)a. (28) 

This result has been obtained on a quite sizable 512 x 512-lattice with the 
help of the cluster algorithm and represents one of the most reliable mass-gap 
correlation length measurements at large values of the correlation length. It 
is consistent with another work In Fig. 9) we display the integration of 
the mass-gap in units of A[ att from its start value at — 1.8 up to the largest 
considered /3-values, the solid curve in the figure. We obtain a determina- 
tion of the mass-gap of the theory in units of the three-loop lattice cut-off 
parameter 

m = 79.62 ± 1.92 A latt . (29) 

The quoted error is of statistical nature. This result agrees with the analytical 
result m = 80.0864A; at4 derived from the thermodynamic Bethe Ansatz. 

5 Conclusion 

The calculation of the mass-gap in the asymptotically free D=2 nonlinear 
sigma model from the numerical evaluation of the path integral poses a non- 
trivial problem. In the standard formulation of the theory asymptotic scaling 
is only exhibited for values of the nearest neighbor coupling /3 > 2.5. Direct 
correlation function measurements of the mass-gap are ruled out. A theo- 
retical device is needed in order to connect the small /3-region of the theory 
with the perturbative regime. 

In this paper we argue, that the consideration of twisted spin configu- 
rations in union with a finite size scaling analysis is able to bridge the gap 
inbetween the small and large /3-regions of the theory. It is possible to de- 
termine the coupling parameter flow of the spin stiffness correlation length 
in a /9-interval which connects both regions. Focussing on the theories spin 
stiffness allows us in addition to adjust certain parameters of the finite size 
scaling law to their perturbative values. These are obtained either in the 
numerical simulation or can be given by theoretical arguments. The mea- 
sured asymptotic value of the spin stiffness correlation length serves as an 
integration constant for the intgration of £ s towards lower /3-values, while the 
spin stiffness constant p s adjusts the amplitude of single logarithmic scaling 
terms of the considered scaling law. This facilitates the precise extraction of 
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the A/3(/3)-shift corresponding to the nearest neighbor coupling dependence 
of the stiffness correlation length. Our analysis strongly supports the exis- 
tence of a scaling window starting at f3- values of about 1.8 with an extension 
into the perturbative scaling regime. Within this region the mass-gap and 
stiffness correlation length are presumed to flow in coupling parameter space 
under the rule of a universal A/9(/9)-shift. The integration of the shift in 
union with a start value for the mass-gap confirms the known theoretical 
mass-gap result. During the course of our study we have not found any 
indication in favor of the existence of a phase with vanishing mass-gap in 
the sigma model. It would be interesting to apply the presented method to 
twisted gauge field configurations and calculate the mass-gap of lattice QCD. 

Acknowledgments: The author likes to thank A. Billoire for helpful 
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Tables 



p 




1.8 


.10055(100) 


1.9 


.10233(051) 


2.0 


.10598(053) 


2.1 


.10966(048) 


1. 1 


.11268(044) 


2.3 


.11486(042) 




.lloz4(040) 


2.5 


.11697(038) 


2.6 


.11724(035) 


2.7 


.11719(034) 


2.8 


.11693(034) 


2.9 


.11655(037) 


3.0 


.11612(041) 


3.1 


.11568(046) 



Table 1: The A/3( / 3)-shift corresponding to changes of the coupling f3 to 
P — f3 — A(3(f3) under the decrease of the stiffness correlation length by a 
factorof 2 : £ a 0) = #,(0). 
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Figure Captions 



Figure 1: Flip-autocorrelation time t/ of the Multicanonical Ensemble 
simulation as defined in section 3 for a fixed L = 36 lattice as a function of (3 
(circles), and at fixed (3 = 2.4 as a function of the linear extent L (triangles). 

Figure 2: Graphical representation of a field configuration on a 36 2 
lattice at vanishing mean field $. The fields on paths are mapped onto the 
complex plane and neighbors within a path are connected by lines. In Fig. 
2a) 36 such mappings are plotted within their corresponding unit circles, if 
the paths are chosen to be the main axis of the lattice in the 1-direction. 
Fig. 2b) contains the corresponding mappings, if the paths are chosen in the 
2-direction. A winding angle of 2rc is observed in the latter case. The (3-value 
is = 8. 

Figure 3: The CEP on lattices of size L = 24 (circles), L = 36 (triangles) 
and L = 70 (crosses) at (3 = 1.6 (Fig. 3a)) and at (3 = 2.4 (Fig. 3b)). 

Figure 4: The quantity AU — 2n 2 /3 on a 36 2 lattice as function of the 
coupling (3. The inlay displays the quantity AU without subtraction. The 
solid curve corresponds to the fit as described in section 4 of the paper. For 
reasons of simplicity we omit the drawing of error bars, except for one data 
point. 

Figure 5: The quantity AU as a function of InL for 12 selected values 
of the coupling f3. The various curves corresponds to the fit as described in 
section 4 of the paper. Note that the data at (3 = 1.6 have not been included 
into the fit. In fact small deviations of the data from the fit become visible 
for the larger lattice sizes at this value of (3. For these data the control 
parameter x s is not large enough. 

Figure 6: The quantity C g as a function of (3, the solid curve. The 
dotted curves indicate the error intervals as supported by the analysis. The 
dashed horizontal line represents the classical limit. 

Figure 7: Logarithmic plot of the stiffness correlation length £ s , as ob- 
tained in section 4 of the paper. The fat circle corresponds to the start value 
for C, s in the asymptotic scaling region. The inlay exhibits the /^-derivative 
of ln^ s characterized by the solid curve. Error intervals are indicated by dot- 
ted curves. The dashed curve of the inlay represents the 3-loop perturbative 
result. The inlay also contains a dotted horizontal line at a value of 2ir, the 
constant term in the inverse /^-expansion of the considered quantity. 

Figure 8: The A/3(/3)-shift of the spin stiffness correlation length (solid 
curve) for values of (3 > 1.8. The dashed curve again corresponds to the 
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3-loop perturbative result. Errors intervals are again indicated by dotted 
curves. The crosses correspond to numerical results of Hasenfratz and Nie- 
dermayer, obtained by a real space renormalization group matching proce- 
dure. 

Figure 9: The integration of the mass-gap in units of the lattice cut-off 
for values of f3 > 1.8, the solid curve in the plot. Error intervals are treated 
as before. The dashed hroizontal line corresponds the exact result. The start 
value for the integration is indicated by the solid square. The asymptotic 
value of the mass- gap is indicated by the solid circle. 
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